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1 Introduction 



It is well known that colliding nuclei break into several small and medium size pieces and 
a lot of nucleons are also emitted. This subfield, known as multifragmentation, gained 
momentum after several theoretical and experimental groups around the world put their 
collective efforts to understand this process [1]. Because of accumulation of experimental 
data on multifragmentation, one has the opportunity to study the role of dynamical 
correlations in fragment formation. Theoretically, the availability of large number of 
models makes the situation worse [1-4] . Due to fact that fragmentation need fluctuations 
and correlations, the molecular dynamical models are the only resource in theoretical 
domains. The molecular dynamics (n-body) approach is well suited as it incorporates 
the correlations and fluctuations among the nucleons. We will, therefore, use quantum 
molecular dynamics (QMD) model [2, 3] to study the dynamics of heavy-ion colhsions. 
Since every model simulate single nucleon, one needs to have afterburner to clusterize the 
phase space. In a very simple picture, wc can define a cluster by using space correlations. 
This method is known as minimum spanning tree (MST) method [5]. In this method, we 
allow nucleons to form a cluster if their centroids are less than 4 fm. This method works 
fine when the system is very dilute. At the same time fragments formed in MST method 
will be highly unstable (especially in central collisions) as there the two nucleons may not 
be well formed and therefore can be unstable that will decay after a while. In order to fllter 
out such unstable fragments, we impose another cut in terms of relative momentum of 
nucleons. This method, dubbed as minimum spanning tree with momentum cut (MSTP) 
method was discussed by Puri et al. [6]. Unfortunately this study was restricted to heavier 
systems like ^^Nb+^^Nb and ^^^Au+^^''Au reactions. The role of momentum cut on the 
fragment structure of lighter systems is still unclear. We aim to address this in present 
paper. 

Exclusively we plan to 

(i) see the role of momentum cut on the fragment structure of lighter colliding systems. 

(ii) and see the role of colliding geometry on the fragment structure with momentum cut 
being imposed. 

The present study is carried out within the framework of QMD model [2, 3] which is 
described in the following section. 
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2 The Formalism 



2.1 Quantum Molecular dynamics (QMD) model 

We describe the time evolution of a heavy-ion reaction within the framework of Quantum 
Molecular Dynamics (QMD) model [2, 3] which is based on a molecular dynamics picture. 
This model has been successful in explaining collective flow [7] , elliptic flow [8] , multifrag- 
mentation [9] as well as dense and hot matter [10]. Here each nucleon is represented by a 
coherent state of the form 

3 

Thus, the wave function has two time dependent parameters and Pa- The total n-body 
wave function is assumed to be a direct product of coherent states: 

= (f)aixi,Xa,Pa,t)(j)l3{x2,Xfi,P)S,t)...., (2) 

where antisymmetrization is neglected. One should, however, keep in the mind that the 
Pauli principle, which is very important at low incident energies, has been taken into 
account. The initial values of the parameters are chosen in a way that the ensemble 
(Arp+Ap) nucleons give a proper density distribution as well as a proper momentum 
distribution of the projectile and target nuclei. The time evolution of the system is 
calculated using the generalized variational principle. We start out from the action 

r c[<t^,r]dT, (3) 



tl 



with the Lagrange functional 

C 



dt 



(4) 



where the total time derivative includes the derivatives with respect to the parameters. 
The time evolution is obtained by the requirement that the action is stationary under the 
allowed variation of the wave function 

5S^5 [' /:[(f),(j)*]dt^O. (5) 

If the true solution of the Schrodinger equation is contained in the restricted set of 
wave function (f)^ (^^i^^a^Pa 

) , this variation of the action gives the exact solution of the 
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Schrodinger equation. If the parameter space is too restricted, we obtain that wave func- 
tion in the restricted parameter space which comes close to the solution of the Schrodinger 
equation. Performing the variation with the test wave function (2), we obtain for each 
parameter A an Euler-Lagrange equation; 

ddC dC 

For each coherent state and a Hamiltonian of the form, 

H = Ta + ^ Kt/sj , the Lagrangian and the Euler-Lagrange function can be 
easily calculated [?] 



^ 2Lm' 



m ^ 

Pa = -Vx„ (9) 

Thus, the variational approach has reduced the n-body Schrodinger equation to a set of 6n- 
different equations for the parameters which can be solved numerically. If one inspects the 
formalism carefully, one finds that the interaction potential which is actually the Briickner 
G-matrix can be divided into two parts: (i) a real part and (ii) an imaginary part. The 
real part of the potential acts like a potential whereas imaginary part is proportional to 
the cross section. 

In the present model, interaction potential comprises of the following terms: 

K/3 = + + Vcoul + VYuk (10) 

Vioc is the Skyrme force whereas Vcoui, Vvuk and Vmdi define, respectively, the 
Coulomb, and Yukawa potentials. The Yukawa term separates the surface which also 
plays the role in low energy processes hke fusion and cluster radioactivity [12]. The 
expectation value of these potentials is calculated as 

^llc = J fa{Pa,ra,t)ff}{pp,r0,t)v}-'^\ra,rp) 

xd'^r^c/'^r^ii^PaC^^P/?, (11) 
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Vilc = J fa{Pa,ra,t)fi3{pis,rfs,t)f^{p^,r^,t) 

xd^Pad^Pfsd^p^. (12) 

where /a(Pa, i"q;, t) is the Wigner density which corresponds to the wave functions (eq. 2). 
If we deal with the local Skyrme force only, we get 



A-r+Ap 
y Skyrme ^ ^ 

a=l 



(13) 



Here A, B and C are the Skyrme parameters which are defined according to the ground 

state properties of a nucleus. Different values of C lead to different equations of state. A 
larger value of C (= 380 MeV) is often dubbed as stiff equation of state. The finite range 
Yukawa (VYuk) and effective Coulomb potential (Vcoui) read as: 

exp{-\ri - rj|}//x 



_ ^ ea;p|-|ri-rj||//x 



Vcoui = V (15) 

The Yukawa interaction (with -6.66 MeV and /i — 1.5 fm) is essential for the 
surface effects. The relativistic effect does not play role in low incident energy of present 
interest [11]. 

The phase space of nucleons is stored at several time steps. The QMD model does 
not give any information about the fragments observed at the final stage of the reac- 
tion. In order to construct the fragments, one needs clusterization algorithms. We shall 
concentrate here on the MST and MSTP methods. 

According to MST method [5], two nucleons are allowed to share the same fragment 
if their centroids are closer than a distance rmin, 

|ri - rjl < rmin- (16) 

where ri and rj are the spatial positions of both nucleons and r^j„ taken to be 4fm. 

For MSTP method,we impose a additional cut in the momentum space, i.e., we allow 
only those nucleons to form a fragment which in addition to equation(16) also satisfy 

|Pi -Pj I <Pmm, (17) 

where p^j„ = 150 MeV/c. 
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Figure 1: The time evolution of A"*"^, free nucleons and LCPs for the reaction of ^^C+^^C 
at incident energy of 100 (left panels) and 400 MeV/nucleon (right) with MST and MSTP 
methods, respectively. 

3 Results and Discussion 

We simulated the reactions of ^^C+i^C , ^oCa+^^Ca, ^SZr+^^Zr and ^^^Au+^^^^Au at 100 
and 400 MeV/nucleon at 6 = 0.0, 0.2, 0.4, 0.6 and 0.8. We use a soft equation of state 
with standard energy-dependent Cugon cross section. 
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Figure 2: Same as Fig. 1 but for the reaction of ^°Ca+^°Ca. 

In Figure 1, we display the time evolution of A'"'^^[(a),(b)], free nucleons [(c), (d)] 
and LCPs(2<A<4) [(e), (f)] for the reactions of ^^c+i^c at 100 (left panels) and 400 
(right) MeV/nucleon. Sohd lines indicate the results of MST method whereas dashed hues 
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represent the results of MSTP method. The heaviest fragment A™"^ follows different time 
evolution in MSTP as compared to MST method. In MST, we have a single big fragment 
whereas momentum cut gives two distinct fragments which shows realistic picture. 

In Figure 1(c) and 1(d), we display the time evolution of free nucleons. We see that for 
both the energies, MSTP method yields more free nucleons compared to MST method. 
There is also a delayed emission of nucleons in MST because of no restrictions being 
imposed. This delayed emission of free nucleons in MST method taken place because of 
the fact that till 30 fm/c, we have a single big fragment in MST method (See Figure 
l(a),(b)). The fragments saturate earher in MSTP than MST as predicted in Ref. [6]. 

In figure 1 [(e), (f)], we display the time evolution of LCPs. We see that MST 
yields more LCPs. The difference between MST and MSTP method increases at 400 
Mev/nucleon signifying significant role of momentum correlations at higher incident en- 
ergies. 

In figure 2, we display the time evolution of A"*"^, free nucleons, LCPs and IMFs for 
the reactions of ^°Ca-|-^°Ca respectively at 100 (left panel) and 400 (right) MeV/nucleon. 
We see that A"*"^ and free nucleons follow similar behavior as reported for the reactions 
of ^^C+^^C. The emission of free nucleons enhances with the cut and thus reducing the 
number of LCPs and IMFs. Similar effects also seen for the reaction of ^^Zr+^^Zr. 

In figures 3 and 4, we display the impact parameter dependence of A™"^, free nucleons, 
LCPs, and IMFs for the reaction of ^°Ca+^°Ca and ^^Zr+^^Zr, respectively, at 100(left 
panel) and 400 (right) MeV/nucleon. From both figures we see that A™"^ rises with impact 
parameter for both methods uniformly. The difference increases with impact parameter. 
This happens because of the fact that we have a bigger spectator matter (from where 
pmax generates) at peripheral coUisions geometry. The number of free nucleons decreases 
with increase in impact parameter for both methods. 

The emission of fragments shows that in contrast to central collisions, peripheral col- 
lisions does not show drastic changes with method. This happens due to fact that with 
increase in colliding geometry, the fragments are the rcmiuauts of either projectile or tar- 
get, therefore breaking mechanisms are almost bound and therefore MSTP methods does 
not give different results. 

In figure 5, we display the impact parameter dependence of A"*"^, free nucleons, and 
LCPs for the reaction of ^^C+i^C at 100 (left panel) and 400 (right) MeV/nucleon. We 
see that in this particular case, the effect of momentum cut on the fragment production 
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enhances with impact parameter (see figure 6(e) and (f)) which is quite different compared 
to earher figures. This is because the spectator matter even at peripheral geometries will 
be very less in such a lighter system and so the fragments are emitted mostly from the 
participant region, where they are unstable and hence momentum cut plays a role at such 
geometries. 

In figure 6, we display the percentage difference [{AX{%) — (X^^mst) — X{rnstp)l Xmst) x 
100)] . We display the system size dependence of the percentage difference of the quantities 
pmax £j.gg nucleons at h— 0.0, 0.2, 0.4, 0.6 and 0.8. We see that percentage difference 
of A"*"^ (circles) is almost constant and remains independent of system size at central and 
semicentral colliding geometry where it is more for lighter systems at peripheral colliding 
geometries. Similar behavior is also observed for the free nucleons. 

In figure 7, we display the system size dependence of the percentage difference of 
LCPs and IMFs at various impact parameters. From figure we see that in central col- 
hsion, AIMF% is almost independent of the system size whereas at peripheral colliding 
geometries it increases with system mass. The average difference is 22, 69, -45, and -16 
for A"*"^, free nucleons, LCPs, and IMFs, respectively. 

4 Summary 

Using quantum molecular dynamic model,we studied the role of momentum correlations 
in fragmentation. This was achieved by imposing cut in momentum space during the 
process of clusterization. We find that this cut yields significant difference in the multi- 
fragmentations of system at all colliding geometries. 
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Figure 3: The impact parameter dependence of A™"^, free nucleons, LCPs and IMFs for 
the reaction of ^°Ca+^°Ca at 100 (left panels) and 400 (right) MeV/nucleon with MST 
and MSTP methods. 
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Figure 4: Same as Fig. 3 but for the reaction of ^^Zr+^^Zr. 
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Figure 5: Same as Fig. 3 but for the reaction of ^^C+^^C. 
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Figure 6: The system size dependence of the percentage difference of A™"^ and free nucle- 
ons between MST and MSTP at various impact parameters. Lines represent the average 
values. 
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Figure 7: Same as Fig. 6, but for the LCPs and IMFs. 
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